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ARTICLE VII. 



New Formulae, relative to Comets. By E. Nutty, Philadelphia. Read 
September 21, 1838. 

The investigation which I here propose to make, respects the com- 
ponent velocities of a comet, observed in three positions, at consecutive 
and moderately small intervals of time. It has for its basis the theo- 
rem of Maclaurin, as adapted to proximate states of a variable function, 
and the knowu expressions for the sun's attractive force on the comet 
and the earth, referred as usual to rectangular solar axes. The means 
which I employ are therefore the same as those presented by La- 
grange in the Mecanique Analytique, and which Mr Pontecoulant 
has recently adopted in his Theorie Analytique du Systeme du Monde, 
where formula? for determining the distances and orbits of comets are 
given with appropriate developments. But the object which I have 
here in view, is not the same as that of Lagrange, in his celebrated 
work above mentioned ; and my investigation and results are different 
from those of Mr Pontecoulant, and embrace a wider extent of sub- 
ject. Similar diversity and extension, in mathematical research, are 
in perpetual requisition. They constitute an essential and important 
part of analytic science ; and with their peculiar attractions, always 
lead to useful views and advantageous contrast. As to the instance 
now adduced, the presumed novelty, and the great accuracy and sim- 
plicity of the formulae which I have obtained, entitle them, I should 
vi. — 3 T 
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hope, to attention and preference. Their mode of investigation I also 
judge important, and as peculiarly eligible. It has enabled me to ex- 
hibit the formulae hitherto given, as particular states of those just no- 
ticed ; and besides others equally simple, it has furnished two new and 
general sets of expressions for the exceptive cases in which the observed 
latitudes and longitudes of the comet would render the general formulae 
doubtful or indeterminate. To the analysis of the principal of these 
results, and with regard to practical applications, I have adjoined the 
data of the comet of 1805, and for which I am indebted to the excel- 
lent treatise of Mr Pontecoulant. The corresponding velocities I have 
computed by the formulae now given, and by others connected with 
the method of La Place. Their comparison has led me to some re- 
marks with which I conclude this paper, and which I have inserted 
from an opinion of their analytical and practical importance. 



Before I enter on the proposed investigation, I think it may not be 
improper to observe, that within a few days, I have been favoured with 
the perusal of Mr Encke's Astronomical, Annual Register {Astrono- 
mische Yahrbuch) for 1833, in which its distinguished author has 
given a full and neat analysis of Dr Olbers' method of determining the 
orbits of comets. The greater part of that analysis I had in fact the 
earlier pleasure of reading in Dr Bowditch's Appendix to the Third 
Volume of his Translation of the Mecanique Celeste. But I had not 
been previously apprized of Mr Encke's remarks on methods which 
differ from that of Dr Olbers ; and in this paper, I would be under- 
stood as having no wish to aim at lessening the predilection with which 
I am now acquainted, and which may be well and reasonably founded. 
My own mathematical partialities, I am not unwilling to avow. I 
entertain them on methods and processes of computation which furnish 
symmetrical and direct results ; and in no slight degree am I favourable 
to that method which is connected with the formulae now occupying 
my attention, and which I have endeavoured to present in such form 
as to merit the approval of Mr Encke, and the author of the Theorie 
Analytique du Systeme du Monde. 

Consider a comet at any point C in its orbit, and let its place at the 
distance r from the centre S of the sun, be determined by the rec- 
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tangular co-ordinates #, y, z. Refer the earth's centre E, to the axes 
of x, y, supposed to be in the plane of the ecliptic, and let (X, Y) de- 
termine its place at the distance R from S. The position of the comet 
relatively to E, will then depend on the values of a; — X, y — Y, z; 
and if we denote its co-ordinates measured from E by pa, p/3, py, we 
shall have in its position C 

x—X-\-?a, y = Y-t-(>(3, z = py (1) 

Accent the different letters in these expressions, in order that they 
may correspond to two different positions of the comet at C, C" ; the 
first being supposed to precede, the second to follow C, at the compa- 
ratively small intervals of time /', t". The co-ordinates C, C", in the 
direction of the axes of x, and at the end of these intervals, will then be, 

x' = X' + p'a'; x" = X" -h P "a" ; (2) 

and corresponding expressions will result in the directions of the other 
axes of y and z. 

The determination of these co-ordinates in terms relative to the in- 
termediate position C of the comet, and to the corresponding place E 
of the earth, may be effected by M'Laurin's theorem, and the known 
differentials 



x 



&H — ~~~ ~~~. j -&-II — 



X 



r»» " ~~ R 



? 



which express the sun's attractive force on the comet and earth. By 
means of that theorem, we have the expressions 

x'=x-x,f-[-ix,r-ixJ' 3 ,Sic., X'=X-X/-HX„r-4X,,/\&c. 

which, in virtue of the preceding differentials, take the usual form 

x' = u'x — v'x„ X'=V'X — Y'X > , (3) 

the assumed coefficients of x, x„ and of X, X,, having in terms of the 
interval /', the following values : 
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i'» r /" t' 3 r f* 

u'=l-* Lr + £ LL-, &c; »' = *'_£ ^ +4 ^-, &c. 

S^ ( 4 ) 

The forms for x", X" in (2) are similar to these, observing to change 
the single into double accents, as respects w, U ; v, V ; and also f into 
— /" ; and like expressions will evidently apply in the directions of the 
axes of y, z. 

Substitute (3) and the similar forms for x\ X" in (2), and then eli- 
minate x by the first of equations (1). There will result the two ex- 
pressions: 

V'(x,—X t ) = Vfpa— p'a' + o'f, V" (x,—X,) = pV'_w"pa + vV, 
in which, for brevity, we have assumed ; 

v/g = (V— v')X,— (IT— w')X, vT = (V''—v")X,+(W'—u")X (5) 

v" v' 
and which, being respectively multiplied by —; , — and added, will 

give 

( x ,-X,)(v'+v")=(^v "- * v' >a-^> + £, 9 'a'+v-eW 

(y,-Y,)(« / +« // )=(^ •'- £ »' ) P«~^ ^ + jp PT'+W+W ^ (6) 

z,(t,'+0=(^"-^t,')pa_^py + H; pV; 

the two last being formed by analogy from the first. 

Expressions apparently more simple than these, might have been 
found from the preceding values ; but the present forms are more con- 
venient for the determination of the component velocities x„y„z„ which 
will become known by means of the earth's velocities X„ Y„ when 
we have expressed the geocentric distances p, p', p*, in terms given by 
observation of the comet in the corresponding positions C, C, C". 
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For this purpose, equate the two values of x, — X /5 which result 
from the expressions preceding (5) ; and from these and analogy we 
obtain 



(?+?)p>--| P v-i i py=o i 



(7) 



in which we have put £ = if" — I', q = j?" — 97'. 

Resolve these equations relatively to p, p', p", with their coefficients; 
and in order to exhibit the results in brief terms, let the coefficients 
which are found to affect £, n be denoted by 



A = 0V— 0% A' = P'y— Py", A" = Py'—(3'y -, 
B = a'Y—a'y", B' = ay"— a"y, B' = a'y—ay'. 



(8) 



These are known from observation of this comet at C, C, C". They 
give the conditional equations 



A/? + A/3' 4- A"/3" = 0, Ay -f- A'y' + A"y" = 0, " 
Ba + BV + B V = 0, By 4- B'/ -+- B V = ; ' 



(9) 



and enable us to express the values now under consideration thus : 

V^W F D ' » /p D ' i>" p ~~ D » ^ ^ 

the common denominator being either of the forms 

D = Aa + AW + AW = B/3 4- B'/3' 4. B"0" ; (11) 

which and (9) will be of immediate use in the simplification of for- 
mulae (6). 

vi. — 3 u 
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Values analogous to these, though very different in their form and 
object, are made the basis of Lagrange's method of determining the 
orbits of Comets. They may be seen in his Mecanique Analytique. 

In consequence of the supposed proximity of the extreme positions 
C, C", to the mean position C of the comet, it is easily conceived that 
the denominator D of (10) is a small quantity, liable to be affected by 
unavoidable errors of observations. We shall therefore eliminate it, 
and employ, instead of (10), the following values: 






(12) 



in which we have assumed 

N = A£ + B>7, N' = A'| + B'*, N" = A"£ + BV (13) 

These quantities are the numerators of (10) ; they give, in virtue of the 
conditional equations (9) and (11) 

N« 4- NV + N V = D£, N0 + N'/3' + N'/8" = D* (14) 

Let now (12) be substituted in the general expressions (6), and let 
us attend to the preceding forms for D£, Dj?. We shall then obtain 
the following values : 



vv+v"e 



, = Y, + £ gv? - J N V + «*) +**£* 



(15) 



w u 

in which we have put w (v 1 -^-v") = —, v" % v. 

v J v v 
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These formulae are well adapted to the proposed determination of 
the comet's component velocities. The generality which we have in- 
tentionally observed in presenting them, has imposed on their last terms 
a complex appearance. But this will immediately vanish, when they 
are limited to the particular state in which approximate results are the 
objects of computation. 

For the purpose here intimated, let us employ the expressions (4) 
and their similar forms with double accents. If the intervals /', t" be 
regarded as small quantities, so that the terms which involve the dif- 
ferentials r„ R, of the distances r, R, may be rejected, we shall have 
for the coefficients of N', N" in (15), the values 

v' ~ f * 3r 3 ) ' v"~ t"^ 3rV' 

in which we shall assume for brevity 

fi a fin 

and to a near approximation r may be supposed equal to unity. 

The factor D in the fourth terms of the components x„ y n we have 
already noticed as being a small quantity ; and we shall presently see 
that £, ri depend on factors which never exceed the radius R. We 

may therefore to great accuracy take w = , „ ■ (1 — ~^t\ tne 

second factor of which, analogously to (16), we shall denote by (i. 
Substitute the expressions (4) just employed in the assumed 

quantities (5) ; and put the factor — ^ = k. We shall obtain, 

in consequence of the minuteness of k and the intervals t', t", the two 
values 

whence there results £" — £', or 

^ = ife(f + 0(X-^X / ); 
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and by substitution, the last term of the component velocity x, becomes 

v'£"-hv"? 



v' + v" 



-s fC* I A.,, 



Values similar to these will evidently result for yj and the last term 
of y, ', and if, for the sake of brevity, we put 

£=X-j(f-*')*i n,=Y-l(t<-t'>)Y,; (17) 

so that the quantities £, v\ may take the form 

Z=:ik(t> + Oe, * = P(*'-M>,; (18) 

we shall immediately perceive that after substitution in formulae (15), 
the factors of £„ y, will disappear ; and that all the terms but the last of 
x„ y„ will be independent of the small quantity £ fe. To eliminate 
£ k from these terms, we must have recourse to the first of the equa- 
tions (10), which in virtue of the values (4) becomes 

t'-ht" . t't\ A£ + B>7 



at" 



k 3r 3 J ~ D 



From this and the preceding values of £, *\ we get \\tt't" = 
ft" Do 

(1 — gpr) • aF+B - ' S ° that if f ° r conciseness we P ut the fac " 



/' t" 
tor 1 — -rr-r = v, and 
3r 3 



the general expressions (15) will take the form 

*=*+&$; mw - v u " a " - ^ (tf) + ^° x ' | . 



K 19 ) 
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the several terms of which are completely determined as respects the 
intervals I', i". It remains therefore only to assign values to the earth's 
velocities X„ Y,; and to the equations (1). 

Denote by L the longitude of the earth when the comet is observed 
in its mean position C. Let is be the longitude of the perihelion of its 
elliptic orbit, and e = sin e the eccentricity corresponding to the mean 
distance 1. The value of «s in 1801 was 99° 30' 5", annual in. = 
-h 1' 2". At the same epoch, e was '01685301, sec. var. = 
— *000041 809. The co-ordinates of the earth in terms of the radius 
vector R and longitude L, are 

X = R cos L, Y = R sin L ; 

and the known expressions for R and the elementary area .described in 
the instant dt, are 



cose 2 Ra dL 

H-esin(L— w)' dt 



RUU3 E' t,,, I* MJ 

= ^ ; — 7t ~ -v 5 R ~TT = COS E. 



If the three first of these expressions be differentiated relatively to 

¥ T 

the time / of which R and L are functions ; and -jr be eliminated by 
means of the last expression ; there will result 

X, = p cos L — q sin L, ) 
Y, = p sin L H- q cosL ; ) 

in which we have assumed for brevity 

COS 8 

p = tanesin (L — «), q = -^-; (20) 

and which, in conjunction with the values of X, Y above given, will 
make known the component velocities and position of the earth. 

The forms which we have here adopted, enable us to express the 
values of £„ *i, in simple terms. We may assume for the coefficients 
of cos L, sin L in these values j 
vi. — 3 v 
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acosa=R — £(*'— t")p, 8sma = ^(f — t")q', 
from which there will result 

3R — (r — t)p 3 sin a v J 

and we shall have instead of (17), the values 

£, = s cos(L — a), »?, = s sin (L — a) ; 

in which a is evidently a small arc and s nearly equal to R. 

A farther simplification of the expressions last given, may be effect- 
ed without diminishing their accuracy, or rendering more complex the 
velocities X„ Y,. Conceive the axis of x, the position of which is 
arbitrary, to be directed so as to form an angle equal to a with the ra- 
dius vector R. From this position which we suppose to be less ad- 
vanced than R, the angle L -r a will then take its origin, and we shall 
have 

£,=«, »7, = 0j X==Rcosa, Y = Rsina; (22) 

the corresponding velocities being 

X, = p cos a — q sin a, 1 

XT • ( ( 23 ) 

Y, = p sin a -+• q cos a, ) y 

in which the first term of Y, is extremely small. 

With respect to the angular quantities in the equations (1), now to 
be considered ; let I denote the geocentric longitude of the comet in its 
position C; and A, the corresponding geocentric latitude. If we also 
represent by p the curtate distance of the comet's projected place in 
the plane of the ecliptic, we shall have, in the present position of the 
axis of £, the following values : 

a s= cos (a+J— L), /3 = sin(a4-J — L), y = tana, (24) 

and the co-ordinates (1), by adding their squares, will give the equa- 
tion: 
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r 2 = R 2 -+- 2Rp cos (/— L) -+- p 2 sec V 



(25) 



which establishes a relation between the curtate distance p and the co- 
met's distance r from the centre of the sun. 

Expressions similar to these evidently apply to the extreme posi- 
tions C, C" of the comet; and in (19), (25) we may, if requisite, 
write p cos "k instead of p, which will then denote the ray drawn from 
the centre E of the earth to the comet at C. 

When the distance r is given, we shall know p by means of (25) ; 
and in such case all the quantities in (19) being determined, we may 
compute the values of the required velocities x n y„ z,. But should r 
be also unknown as well as p; and this is generally the case, the equa- 
tion (25) will not alone be sufficient, and we must join to it the diffe- 
rential expression 



x*-\-y*-\-z* = * , 



(26) 



which we adapt as a known form for the square of the comet's veloci- 
ty in a parabolic orbit; and which will enable us to obtain the dis- 
tances r, p, and lastly the component velocities x„ y„ z,. 

We may now present the expressions (19) in their final terms for 
computation, They become, in terms of the simplified values of |„ 
J7„ of the following form : 



--= x -+x(F A ' a '-7 AV -" D T^ 



ox,, 

3 s J' 



y (A) 



in which the values of A, A', A" from (8) are 

A = /?y — (3'Y, A' = /3V — /V, A" = ^y'-/3V; 

and the corresponding value of D is by (11), 

D = Aa-f-AV + AV. 
3 v* 
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There is but one remark to be made relatively to the numerical con- 
version of these and the auxiliary quantities fi, p', &c. and a, s. The 
intervals f, t" and also f — t", given in days, must be expressed in 
terms of mean solar time, and in parts of radius, by multiplying them 
by the diurnal factor -01720213, which corresponds to the solar arc 
59' 8"-2, and of which the logarithm is 8-2355821. 

The preceding formulae (A) are we believe the most convenient 
and the most accurate that have appeared for the determination of the 
component velocities of a comet, observed in three positions at compa- 
ratively small intervals. The factors ^ ft', &c. have here been first 
noticed ; they give every requisite degree of precision to the values of 
x„ y„ z,. The computation we have reduced to uniformity and faci- 
lity by means of the three arcs a-\-l — L, a-\-l' — L, a-\-l"— L, 
of which the several angular quantities a, a', &c. are functions ; and 
which take place of nine arcs hitherto employed. The particular form 
given to D is also advantageous. It directly leads to the value of this 
factor, by the three previously determined quantities A, A', A". 

When the observations can be taken so that the intervals f, t" may 

t 12 
be equal, the factors [i',[i", v will each be expressed by 1 - 3 . The 

term which depends on t'—t" will then vanish, and we shall have 
the more simple expressions : 

v v _i_ p f*' ( A'«' A" " _i- DX ^ 

y ' = Y ' + A? (A ' /3 ' ~ A "^" + ll)' \ (A ') 

*, = j^(Ay-AV); 

in which the quantities A, A', A" and D have the same forms as be- 
fore ; but instead of (24), we have 

a = cos (I — L), = sin (/ — L), y = tan X, 

and similar values for a', a", &c. 
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In the formulae hitherto given for this particular case, the factor yi 
has not been introduced or noticed; and on this account, and the last 
terms involving the earth's velocities X„ Y„ they are less exact than 

DY 

the preceding expressions (A'), in which the term —=£■ at least should 

be computed as no negligible, minute quantity. 

The general formulae, of which we have now completed the inves- 
tigation, and also (A'), fail in giving accurate results, when the quan- 
tities A, A', A" and D are minute ; and they become indeterminate, 
when these quantities vanish. If we conceive the apparent path of 
the comet during the interval /' -|- t" to be a great circle, passing 
through a point in the ecliptic of which the longitude is less than that 
of the sun by the small arc a, the tangent of the inclination of the orbit 

y iy' y" 

will be evidently expressed by either of the three relations ^-, —^ — tt \ 

and we shall have then A = 0, A' = 0, A" = 0, and also D = 0, to 
which the formulae (A) and (A') are inapplicable. The case in which 
these quantities become small, is now obvious without further reflec- 
tion ; and should it be inconvenient or impossible to take new observa- 
tions of the comet at other places in its orbit, we must abandon for- 
mulae (A), and have recourse to different expressions, capable of deter- 
mining the component velocities x„y„ z,. 

Such expressions are the following : they have not to our knowledge 
been hitherto presented. 

, = Y, + | (£ W- |W) + K(f- + ™*) fcB) 

*-6(p»V-t'"V)+K^. 

There are in fact two sets of formulae, comprehended in these ex- 
pressions. We may take in terms of the geocentric longitudes ; 
vi. — 3 w 
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E== sin (/"-/'), E' = sin (/-/"), E" sin (/'-/), e " = (3"; 

or when these longitudes vary from each other less than the corre- 
sponding latitudes, we may more suitably adopt the values 

E = a'y" — a'Y, E' = a"y - ay", E" = ay 1 — ay 1 j e" = y", 

the three first of which are the values of B, B', B" (8) taken with a 
change of signs. In both cases here implied, the values of a, a', &c. 
will be as in formulae (A); and we have taken = t-ht", 

The formulae (E) are regular in the composition of their terms, and 
but little more complex than (A), from which they chiefly differ in 
this respect, that the factor k cannot be eliminated, but must corre- 
spond to a determinate value of the comet's distance r. Any change 
in the value of this distance will however be attended only with small 
additional computation; since the calculated values of the different 
terms in (E) will, from the form which we have adopted, continue in- 
variable. 

When the intervals t, t" happen to be equal, we shall have K = 
Ms?, and the more simple values 

« = x, + £ crv _ e v) + k (■ v _ 4 + vy 

* = T. + qj; m- E"/J«) + K C-| H- 1£) )■ (E) 

the quantities E, E', E", e" being as in (E), and the values of a,a,&c. 
the same as in the particular forms (A'). 

The investigation of the formulae now given, we reserve for another 
paper, which will soon appear. At present we shall briefly show how 
the expressions given by Mr Pontecoulant may be readily deduced 
from our expressions (19). 

If in these we omit the last terms depending on the earth's veloci- 
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ties X„ Y, ; and put p'= 1, (i"= 1, &c. and also p = AMN", there 
will result for the component velocities ar„ y„ z„ the less approxi- 
mate values 

x, = X, + A (/'M'a' — <'MV — (/' — t") D£),l 
y, = Y, + A (* W — *"M"0" - (f — /") D*,), )► (a) 

z, = A (fMV - f My) 5 J 

and if to these be added the identical equations 

Ma + M'a' 4- M"a" — D£ = 0, 
M/S-f-M'/S'-j-M'^" — By,, = 0, 
My + My -j- MV =0, 

respectively multiplied by i h (t" — t'), we shall immediately obtain 
the following expressions : 



a- = X, + h H (MV — M'a") — | & (Ma + D$)}" 

y == Y, 4- h {£ B (M'/3' - M"/3") -KCM^+D.?,)} 
z, = A {£ (My— M"y") — £ & My} ; 



\ 



(*) 



in which = J 7 4*", 0' = /' — /". 

These are in fact less simple than the preceding. They become 
identical with the values F„ G„ H„ found at page 44, vol. II. Theo- 
rie Analytique, &c. when we further neglect the arc a in the values 
of £„ *i„ multiplied by the factor D. The term D£, may be removed by 
directing the axis of y to the earth when the comet is at C. The 
preceding values of x„ y will then agree with P,, Q„ page 45. But 
even in this simplified state, the computation depends on nine different 
angles instead of the three involved in formulae (A). 

The preceding values (a) and (b) are not the only forms that can 
be derived from (19), in virtue of the identical equations which we 
have just employed. If to (a) we add those equations respectively 
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multiplied by h {t" — f ), the terms which depend on D will disappear, 
and we shall obtain a new set of values: 



X, = X, -f- h {t"Ma' — t'M"a" — (/'— /")Ma}, 
y, = Y, -+-h{t"M'P' — t'M'W— (*'— /")M^}, J, 
2, = A {/"My — /'My- (/'— /")My};J 



(0 



which are perfectly symmetrical, and preferable to (J), both as to form 
and facility of computation. 

We presume that the advantage of our mode of solution is suffi- 
ciently tested by the different results we have obtained. We shall 
therefore proceed to the numerical application of formulae (A). In 
the example chosen, the intervals /', /" are considerably different, and 
t — /" is of no small magnitude. We have taken it as before men- 
tioned from the Theorie Analytique du Systeme du Monde, vol. II., 
p. 68. 



Times of Observation. 

Nov.23 d -32241 

" 30-51095 

Dec. 5-29581 



Data of Comet of 1805 

Longitudes Observed. 

I' =24°41 / 04" 



I — 15 39 40 

/" = 2 7 11 



Latitudes Observed. 



W =27°25 / 35" 

a = 19 25 28 

%" = 3 20 45 



Intervals. 



/' = 7 d -18854 
t" = 4 -78486 



These have been corrected for aberration and parallax, and with re- 
ference to mean time at Paris. 

From the tables of Delambre have been taken 

(Long, of + 180°) or L = 68° 25' 41 " ; log R = 9-9936673. 

To these I add from values before noticed ; 

vt == 99° 34' 13"; L — xs = 31° 8' 32"; sin e = -01685151. 
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t' =9 0922228 

«"=5-9154514 

t'— <"=8-6164587 



Preliminary Quantities for (A). 



a + 1'— L=— 42°55'51" 
a+l— L=— 51 57 15 
a+l"—L=— 65 29 44 



p = 


7-9403288 


7 = 


00062710 


a= 


+ 48'46" 



S = 9 9937738 



The numbers in this and the following are sufficiently distinguished 
from the logarithms by the prefixed signs 4- and — 

From Computation of Formulae (A). 
Values ofX, Y,. 



p cos a — q sm a 
•0087163 — -0143914 
X, = — 0231077. 



p sm a + q cos a 

— 0001236+ 10144421 

Y, = + 1-0143185. 





Values of A, A', A", &c. 




p'y" — (3"y' 


p'y-fy" 


fr — py 


— -0398195 + 4722050 


—•3208804 + 0460401 


—•4086707 + -2401851 


A = + 4323855. 


A' = — -2748403. 


A" = — -1684856. 


At' 


At!' At't" 


Aa 


8-7280941 8'< 


5513227 7-64354! 


>5 + 2664756 


/A'a' A"a"\ 
P KJJ 7 ~ At' J 


- /A'/3'_ A"p\ 
*AAf At' J 


/Ay A'y'N 

v \ At" ~ At' ) 


0-7523734; 01162696 


0-7209781; 0-4574862 

1 


06028708; 9 -2653444 


—5-654229 + 1-306982 


+5259907 — 2-867386 


— 4007474 +,0184223 


— p (4-347247). 


+ p (2-392521). 


— p (3-823251). 



We have supposed in this computation that the factors (i, p', &c. are 
each uuity. But if we take for greater precision the distance r = 1 , 

logfi"=9-9990181; log ft'= 9-9977807; log^=99992627; logv = 99985244, 
vi. — 3 x 
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we shall find the corresponding values, 

— p (4-341140). + p (2-395257). — p (3-815140). 

If to the value of A a = -2664756, we join AW = — -2012315, 
A "a" = — -0698817, obtained from the logarithms used in the pre- 
ceding computation, there will result D = — -0046376, and then we 
shall get 



p( D '-ir + D §5) 

8-6392065; 5'9232877. 

+ 0435719 + 0000838 

+ p (-043655). 



(°Hi) 



P 

7-5657013 
— p (-003679). 



In these values we have taken p = 1, fi' == 1, &c. If we employ 
the logarithms which correspond to r = 1, there will result, instead of 
the preceding : 

+ p (043581), — p (003662); 

by virtue of which and the value before given, we obtain for the co- 
met's velocities : 

z, = —0-023108 — p (4-297559), ~\ 

y,= + 1 014319 + p (2-391595), V (d) 

z, = — p (3-815140);, 

in which p is the curtate distance from the centre of the earth. 

Had we retained the values found in case of p == 1, p' = 1, &c, 
the coefficients of p would have been 

(4-303592), (2-388842), (3-823251); 

which are not so accurate as those in (d). 

The preceding expressions for x„ y„ z„ may be easily changed so as 
to correspond to any required position of the axis of x, and to the ray 
drawn from the centre E of the earth to the comet at C. If we sup- 
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posed this axis to be directed to E, we shall find in terms of p before 
used ; 

X, = — 0-008716 — p (4-263202), "j 

y, = + 1-014544 + p (2-452316), \ (e) 

*, = — p (3-815140)., 

But if we also change p into p cos A, so that p may denote the ray 
CE, we shall get 

x, - —0-008716 — p (4-020544), 

y, = + 1-014544 + p (2-312732), 

z, = — p (3-597985). 

The accuracy of these results we believe to be very considerable. 
Mr Pontecoulant's values in terms of p are the following: 

x l = — 0-008686 — p (4026273), y, = + 1-014545 + p (2-314020) ; 

z,— — p (3-605632). (See Theorie Analytique, &c, vol. II., p. 70.) 

With the desire of making a comparison between the principal terms 
of formulae (E), and the values (d) above given, I have subjected them 
to computation. By taking the first expressions for E, E', E" and e", 
the coefficients of the curtate distance p were as here given : 

x, ~ X,—p (4-050920) — &c. 

y, = Y, + p (2047392) + &c 

z, = — p (3-646140) — &c. 

so that the terms affected by K differ little from — p (-245653), 
+ p (-345974), — p (-168997), and will in the destined use of (E) 
not exceed the order of the intervals /', /". 

To the data of the example now considered, I have been induced to 
apply the method of La Place; and with surprise I found results 
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which it may he proper here to insert. By means of the two expres- 



sions 



which are deducihle from the values of l\ 1% expressed by Maclaurin's 
theorem ; I obtained in parts of radius, dl — — 2*232849; d 2 l = 
— 15*51 507. The similar values which depend on the geocentric 
latitudes were d a = — 2-498243, & % h = — 22-139036 ; and 
adopting in this method the formulae : 

i e = -,, »_#| +to (I._j ) (» +iBS __ + _ p ), 
n=2(» + 9tan(L — JW ' U 

1 v y en 



sin a cos ^ 



which may be found expressed differently in the Mecanique Celeste 
and Theorie Analyt 
The three equations 



m 
and Theorie Analytique, &c. I obtained — or * = — 3-278562. 

71 



x, = X y + p (t'a — @dl), 
y l = Y, + p(ip+ ad I), 

resulting from the differentiation of (1), and in which a = cos (/ — L) 
$ = sin (/— L), y = tan a, then gave me the following values : 

x, = X, — p (3-761473), 

y, = Y, + p (1-261091), 

z, = — p (3965141); 

which are exceedingly different from (e), with regard to the coeffi- 
cients of 9 in the values of x, and y,. 
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In addition to these results, I would join the remark, that the com- 
puted values of x„ y„ z„ expressed by the preceding differential for- 
mulae, cannot be made to agree with (e), so long as are employed the 
above numerical quantities found for dl, d'A,; and that from this and 
other instances, I consider as defective the manner in which the first 
and second differentials of the geocentric longitudes and latitudes are 
determined in the method here noticed; and which, I think, should 
not be used for determining the component velocities of a comet, or 
the position of its orbit. 

We might now particularize further the expressions (c?) by means 
of (25) and (26), and the consequent values of the distances r and p. 
But our principal aim in the present paper has, we imagine, been suf- 
ficiently attained. We intend again to resume the subject on an early 
occasion. 
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